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Abstract 

The stabilization of host-symbiont mutualism against the emergence of parasitic individuals is pivotal to the evolution of 
cooperation. One of the most famous symbioses occurs between legumes and their colonizing rhizobia, In which rhizobia 
extract nutrients (or benefits) from legume plants while supplying them with nitrogen resources produced by nitrogen 
fixation (or costs). Natural environments, however, are widely populated by ineffective rhizobia that extract benefits without 
paying costs and thus proliferate more efficiently than nitrogen-fixing cooperators. How and why this mutualism becomes 
stabilized and evolutlonarily persists has been extensively discussed. To better understand the evolutionary dynamics of this 
symbiosis system, we construct a simple model based on the continuous snowdrift game with multiple Interacting players. 
We investigate the model using adaptive dynamics and numerical simulations. We find that symbiotic evolution depends on 
the cost-benefit balance, and that cheaters widely emerge when the cost and benefit are similar in strength. In this 
scenario, the persistence of the symbiotic system Is compatible with the presence of cheaters. This result suggests that the 
symbiotic relationship is robust to the emergence of cheaters, and may explain the prevalence of cheating rhizobia in 
nature. In addition, various stabilizing mechanisms, such as partner fidelity feedback, partner choice, and host sanction, can 
reinforce the symbiotic relationship by affecting the fitness of symblonts in various ways. This result suggests that the 
symbiotic relationship Is cooperatively stabilized by various mechanisms. In addition, mixed nodule populations are thought 
to encourage cheater emergence, but our model predicts that, in certain situations, cheaters can disappear from such 
populations. These findings provide a theoretical basis of the evolutionary dynamics of legume-rhizobia symbioses, which 
Is extendable to other single-host, multiple-colonizer systems. 
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Introduction 

Symbiosis is an ecological interaction in which two or more 
species exchange mutual benefits. This cooperative interaction 
promotes the fitness of both species, and thereby reinforces their 
symbiotic relationship. However, symbiotic systems are vulnerable 
to emerging selfish cheaters that extract benefits fi-om the system 
without paying costs. Such parasitic cheaters may potentially 
disrupt the symbiotic relationship. 

One of the well-known symbioses occurs between legume plants 
and rhizobia (nitrogen-fixing soil bacteria). Rhizobia establish 
symbiotic organs termed root nodules on the roots of their host, 
and proliferate by extracting nutrients from the host plant. In turn, 
they supply their host plants with nitrogen resources produced by 
nitrogen gas fixation. This mutual nutrient exchange should 
promote the fitness of both organisms and thereby strengthen the 
symbiotic relationship. This beneficial effect, known as "partner 
fidelity feedback", is assumed as a stabilizing factor for the 
mutualistic relationship [1-4]. 

On the other hand, naturally occurring rhizobium strains vary 
in their nitrogen fixation activity, and ineffective rhizobia that 
colonize their host plants without undertaking nitrogen fixation in 



their root nodules are ubiquitous [5-8]. Because the nitrogen 
fixation reaction consumes much energy (or costs), such parasitic 
cheaters could use surplus energy for their own growth or for 
synthesizing storage substances. Consequendy, they are likely to 
proliferate more efficiently than nitrogen-fixing cooperators, 
posing a risk to the symbiotic interaction. Rhizobia are therefore 
exposed to two opposite effects that simultaneously promote (by 
providing benefit) and destabilize (by incurring cost) the mutual- 
istic relationship. In this paper, we refer to these effects as the 
"promoting force" and "destabilizing force", respectively. 

Despite the widespread presence of ineffective rhizobia, the 
legume-rhizobia symbiosis is evolutlonarily stable. Potential 
mechanisms that stabilize the symbiotic interaction include 
partner fidelity feedback [1-4], partner choice [4,9-12], host 
sanction [4,10,13-16], spatial structure [17], and kin selection (or 
inclusive fitness) [18]. Partner choice and host sanction exert 
similar effects, in that host plants discriminate beneficial symbi- 
onts. However, partner choice is often regarded as a selection 
process based on recognition signals, while in host sanction, plants 
punish more parasitic cheaters by reducing nutrient supply based 
on their symbiotic performance [3,4,10]. Thus, in this paper, we 
consider partner choice and host sanction as pre-infection and 
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post-infection processes, respectively, according to the nitrogen 
fixation activity of symbionts. 

The evolution of cooperation has been theoretically modeled by 
game theory, in particular, the prisoner's dilemma game and the 
snowdrift game (also known as the hawk-dove game) [19,20]. If the 
game involves two players, mutual cheating is an evolutionarily 
stable strategy (ESS) in the prisoner's dilemma, while the snowdrift 
game has the Nash equilibrium of mixed-strategy, in which 
cooperator and cheater stably coexist. These games have been 
extensively studied under various conditions. However, although 
the simplified framework of binary strategy choice (full coopera- 
tion or fuU defection) captures the essential nature of cooperation, 
evolutionary dynamics involving quantitative traits (such as 
nitrogen fixation activity in rhizobia) require more complex 
frameworks. Continuous strategies, such as the continuous 
versions of snowdrift and prisoner's dilemma, more accurately 
reflect real-world phenomena [20—27]. Interestingly, the contin- 
uous snowdrift game permits an evolutionary process in which 
completely non-productive cheaters coexist with cooperators 
making maximum investment [23]. This evidence provides insight 
into the evolutionary origin of ineffective symbionts in the legume— 
rhizobia mutualism. 

Various mathematical models have been proposed for the 
evolution of the legume-rhizobia symbiosis. Akcay et al. [28,29] 
developed models based on bargaining theory, in which legumes 
and nitrogen-fixing nodule(s) negotiate to d(x ide whether and how 
to form or terminate the interaction. Other models describe the 
population dynamics of rhizobia strains with fixed nitrogen 
fixation activity, such as complete cooperators or total cheaters 
[20,30-33]. These models can explain the coexistence of 
cooperators and cheaters; however, as described above, discrete 
strategies oversimplify the evolutionary dynamics of quantitative 
phenotypes. On the other hand, West et al. [14,18] reported 
models in which the optimal strategy (or the maximum fitness) is 
determined by continuously varying the nitrogen fixation. While 
finding the optimal strategy is a powerful approach for analyzing 
evolutionary systems, it cannot sufficiently explain relatively 
complicated situations, such as frequency-dependent interactions 
between different rhizobia strains. In fact, the stable co-existence 
of cooperators and cheaters is difficult to explain using these 
approaches. 

In contrast, the legume-rhizobia mutualism is more directiy 
reflected by evolutionary dynamics in which rhizobia prohferate 
according to the quantitative trait of nitrogen fixation with a 
frequency-dc'pcndcnt selection. However, this situation has yet to 
be fully investigated, although such a model may explain the effect 
of mixed colonization [34]. Therefore, to theoretically understand 
the evolutionary dynamics of the legume/ rhizobium system, we 
construct a simple model based on the continuous snowdrift game 
with interacting multiple players [20,23], and investigate the 
model using adaptive dynamics and numerical simulations. The 
aim of this paper is not to examine the nature or evolution of 
mechanisms that prevent rhizobia from cheating, but rather to 
investigate the conditions under which mutuaHstic and cheating 
symbionts can coexist. We also aim to understand how and why 
symbiosis establishes and stably persists despite the ubiquitous 
distribution of cheating rhizobia. 

Models 

2.1 Model Framework 

2.1.1 Fitness function. In the legume-rhizobia symbiosis, 
host plants cannot extract benefits (i.e. nitrogen fixation products) 
without paying costs (i.e. nutrients or energy) to rhizobia, because 



the symbionts are confined in root nodules and thereby cannot 
directly interact with the external environment. This fact suggests 
that host plants cannot develop a cheating strategy, whereas 
ineffective rhizobia are ubiquitously distributed in nature. Thus, 
the evolutionary dynamics of the mutualism are most likely driven 
by the strategies of the symbionts, while the host plants play a 
minor role. Accordingly, our model is based on the strategies 
adopted by the rhizobia. 

The evolutionary strategy of each rhizobium is its nitrogen 
fixation activity in a root nodule, which is continuous between 0 
(no activity or full defection) and 1 (maximum activity or full 
cooperation) and which is genetically transmitted to its progeny. 
Each host plant is assumed to be colonized by a constant number n 
of rhizobia, which correspondingly form n root nodules (Figure 1 A). 
Thus, each root nodule harbors a clonal monomorphie population 
of rhizobia, except in the mixed-population model of section 3.8, 
where each nodule harbors one or more than one rhizobial strain. 
In addition, host plants are assumed to be stochastically infected 
with rhizobia independent of their nitrogen fixation activity (or 
strategy), except in the partner choice model of section 3.7, where 
host plants select preferential symbionts as colonists. 

Because rhizobia are considered to grow much more efificientiy 
in root nodules than in the soil, we consider only their proliferation 
in root nodules. Proliferation is linked to the benefit (or promoting 
force) and the cost (or destabilizing force) of nitrogen fixation 
activity. The benefit encourages cooperative growth between the 
host plant and its symbionts. Now consider that a host plant is 
infected with n rhizobia of nitrogen fixation activity (^= 1, 2, 
— , «) according to the probability distribution of rhizobia in the 
soil, and subsequently generates n corresponding root nodules. 
The productivity of this focal plant {(p^ is promoted by the total 
nitrogen resource supplied by the root nodules, represented by the 

n 

average of their strategies: x= ^ Xkjn, and is thus an increasing 

4=1 

function of x: q)j,=B\{x). Each colonizing rhizobium (or its 
corresponding root nodule) in turn grows more efficientiy as the 
host plant productivity (pp increases. Thus rhizobial growth rate (^) 
is an increasing function of (pp. (p = B2{(Pp), and is consequentiy 
expressed as (j> = B2{B\{x)) =(PqB(x) and B(0) = 1, where q>o is the 
growth rate in the absence of nitrogen fixation (x = 0) and the 
benefit function B[x) is an increasing function of x. 

On the other hand, the cost (or destabilizing force) exerts an 
inhibitory effect on symbiont fitness, such that the growth rate of 
each infecting rhizobium (or its corresponding root nodule) 
decreases as its own nitrogen fixation activity (x,) increases; 
^, = ^gC(x,) and C(0)=1, where the cost function C{x) is a 
decreasing function of x. Therefore, when benefit and cost 
compete, the growth rate of a focal infecting rhizobium adopting 
strategy can be expressed as ^, = ^q-6(x)C(x,). The transfor- 
mation of = g>Qq)j yields the relative fitness: 

n 

p,=5(x)C(x,),5(0) = C(0) = l,andx= ^Xk/n, (1) 

k=\ 

where the benefit and cost functions B{x) and C{x) are strictiy 
increasing and decreasing functions ot x, respectively, (i.e. B'(x)>0 

and C'(x)<0). 

Because the promoting and destabilizing forces affect rhizobia 
proliferation in a synergistic rather than an additional manner, we 
define the fitness as the product of the benefit and cost functions 
(Eq. (1)). This form differs from that of conventional game theory, 
in which the cost is subtracted from the benefit. However, 
qualitatively, both forms of the fitness function yield the same 
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Figure 1. Schematic representation of model framework. (A) The evolutionary strategy of an individual rhizobium is its nitrogen fixation 
activity (Osxsl). Rhizobia stochastically inhabit host plants according to Pt(x), the probability distribution of their strategy, proliferate in root nodules 
in a frequency-dependent manner according to Eq. (1), change in their nitrogen fixation ability by mutation, and are released back to the soil 
following the death of their host plants (Eq. (4)). Proliferation of colonized rhizobia is driven by the benefit (promoting force) and cost (destabilizing 
force), which depend on the nitrogen fixation activities of the rhizobia. The cost of a single rhizobium is affected by its own strategy x, (i.e. C(x,-)) while 
the benefit is affected by x: the average strategy of n colonizers (i.e. B{x)). See text for details. (B and C) Effects of partner choice (B) and mixed 
nodule (C) on the fitness of rare mutants with nitrogen fixation activity y in the resident population with x. a(x) is the probability that a rhizobium with 
X is accepted by a host plant, and fi is the frequency of nodules inhabited by two rhizobium species. Yellow and blue root nodules are colonized by 
residents with x and rare mutants with y, respectively. 
doi:1 0.1 371 /journal.pone.0093670.g001 



evolutionary dynamics of the symbiotic system, because the 
evolutionary consequences are similar in the two models; namely, 
full defection, fuU cooperation, partial cooperation, and coexis- 
tence of cooperators and defectors created by "evolutionary 
branching" (corresponding to cases (i), (ii)/(vi), (iii), and (iv)/(v), 
respectively, in this paper) [23]. Consequently our model is 
equivalent to the continuous snowdrift game with n interacting 
players [20,23]. 

2.1.2 Benefit and cost functions. The benefit and cost 
functions are simply defined as: 

B(x) = l + bx- bf^hx^ and C(x) = 1 - cx + c^cx^ , (2) 

where b and bj^ correspond to the strength and nonlinearity, 
respectively, of the benefit, and c and are the corresponding 
parameters of the cost (Figure lA). These functions require that 
B(x)>0, B'{x)>0, C(x)>0, and C'(x)<0, which are satisfied when 
0<A, 0<e<l, 0<A,v<l/2, and 0<Cjv=^l/2. In this paper, C{x) is 
not regarded as an increasing function of x, because nitrogen 
fixation obviously strengthens if both the benefit and cost functions 
are increasing functions of x (see Text S2). 

Partner fidelity feedback and host sanction are widely-discussed 
stabilizers of symbiosis [2-4,35-37], and are incorporated into the 



benefit and cost functions, respectively, in our model. The effect of 
the benefit (/>) increases as partner fidelity feedback strengthens, 
while the cost (c) decreases with increasing host sanction. 

2.2 Numerical Calculations 

In the numerical calculations, we first introduce the discrete 
probability distribution of the rhizobia strategies in the t-th cycle: 

N 

and ^p,(/)=l, (3) 

i = 0 

where 2 (0, 1, — , jV) corresponds to strategy x = i/JV. The initial 
distribution is a population of non-fixing rhizobia: Pq{0) = 1 and 
po({) = 0 for zT^O. This strategy distribution is iteratively changed in 
the infection, proliferation/ mutation, and release steps (Figure 1 A). 
We consider a total number A4 of host plants, each infected with n 
rhizobia randomly selected from the probability distribution pit). 
We denote the probability distribution of the Mxn infecting 
rhizobia by p'/ii). In each corresponding root nodule, the infecting 
rhizobia proliferate at a rate calculated by Eq. (1) and their 
nitrogen fixation ability is altered from strategy J/ jV to ?/jV by 
mutation at transition rate pj,. Thus the probability distribution 
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after the proliferation/mutation step can be described as a time- 
discretized form of the rephcator-mutator equation: 

Pt+ 1 (0 = X] f^ji9>jP't(j) / 9, (4) 

7 = 0 

N 

where (p= ^ VjP'tij) is the average fitness of the infecting 

7 = 0 

rhizobia, and (pj is the fitness of strategy j/N. In this paper, 
mutation rates are given by 

tiji = /io ( = 0) ,/^/, = ( / =7 + 1 ) , 

Hji = 0{i^OJJ± l),and/i,., = 1 - ^Mk, 

where /i specifies the mutation rate between neighboring 

strategies, and /.(o specifies the rate of the null mutation. 
Parameters are set to jV= 200, M=2000, ^-0.04, and fio-O.O 
or 0.01. 

Results 

By investigating our model, we hope to understand the 
evolutionar)- ronsc^quences of the legume-rhizobia symbiosis in 
terms of adaptive dynamics (Text SI) [38,39]. Adaptive dynamics 
assumes that strategy changes induced by mutations are small. 
However, this assumption may not hold in legume-rhizobia 
symbiosis, because a defect in a single gene can impart the null 
mutation that is responsible for complete lack of nitrogen fixation 
activity. Thus, we first analyze the model assuming that mutations 
exert litde effect, and afterwards examine the impact of null 
mutations on the symbiotic relationship. 

In our model, a host legume plant is infected with a constant 
number n of rhizobia, each of which adopts a nitrogen fixation 
activity O^x^l as its evolutionary strategy (Figure lA). Now 
consider the invasibility of mutants with nitrogen fixation activity jf 
in the resident population with x. If the mutant is rare, the 
probability that a host plant is colonized by multiple mutants is 
very low and can be neglected. Thus, assuming that a host plant is 
colonized by a single mutant withy and by M - 1 residents with x, 
Eq. (1) gives the fitness of the mutant strain as 

f,(y) = Bix)C(y)andx = i{n- \)x+y)/n. (6) 

Note that, as described in section 2.1.1, the fitness is the product 
of the benefit and cost functions, unlike conventional game theory, 
in which these quantities are subtracted. We now introduce the 
relative fitness of the mutant to the resident: 

Hx,y) =fM -fxix) = B{x)C{y)-B{x)C{x), (7) 

which determines the invasibility of the mutant into the resident 
population. \i w(xJ]>Q, the mutant can invade the population, but 
cannot if w(xy)<0. The theoretical basis of adaptive dynamics is 
described in Text SI. 

Because legumes and rhizobia provide mutually-enhancing 
nutrients to each other, thereby increasing the fitness of both, the 
benefit function B[x) should be a strictiy increasing function of x 
(i.e. B'{x)>0). In contrast, the cost function C(x) is not necessarily a 
simple decreasing function of x because it may be affected by 



partner choice and host sanction. If such stabilizing effects are 
suflScientiy strong that Qa;) becomes an increasing function of x, 
the maximum activity of nitrogen fixation [x=\) is predicted to 

evolve and be stably maintained (see Text S2). Accordingly, we 
assume that C[x) is a decreasing function of x (i.e. C'(x)<0) in this 
paper. Note that under this condition, more cheating rhizobia 
proliferate more eflSciently when they inhabit the same plant. 

3.1 Linear Cost Function 

In the linear cost function (i.e. c^- 0), the selection gradient 
D{x) = dfx(y)ldy\y,^^ is a decreasing function of A' (i.e. D'{x)<0 for 
0<x<l). Accordingly the evolutionary behavior is classified into 
the following three cases depending on the cost-benefit balance 
(for details see Text S3). 

3.1.1 Case (i) "No evolution". If Z)(0)<0 (i.e. c>Co = */«), 
we have Z)(x)<0 for (i<x<\ (Figures 2A and SI, gray). In this 
scenario, a resident population with nitrogen fixation activity x is 
always invaded by lower-performing mutants (y<x) but not by 
those with higher activities {y>x) (Figure 2B), and a symbiotic 
relationship never evolves. We refer to this scenario as the "No 
evolution" case. 

3.1.2 Case (ii) "Maximum evolution". Conversely, if 
Z>(1)>0 (i.e. c<ci={\-2hM)h/{n + {n+\-nbM-2hM)b}) 
(Figures 2A and SI, magenta), we have Z)(x)>() for ()<jc<l. In 
this scenario, the nitrogen fixation of the resident population 
always increases, because only mutants with higher strategies can 
invade and subsequently replace the resident (Figure 2C). There- 
by, the population evolves toward the maximum nitrogen fixation 
of x= 1. We refer to this scenario as the "Maximum evolution" 
case. A population that has evolved to x = 1 cannot be invaded by 
cheating bacteria (x = 0) possessing the nuU mutation (because 
41,0)<0; see Text S3.4). 

3.1.3 Case (iii) "Intermediate evolution". The interme- 
diate condition between cases (i) and (ii) (i.e. c\<c<q^ (Figures 2A 
and SI, blue), permits a unique singular strategy 0<x <1 that is 
always convergence stable (CS) and ESS-stable (i.e. D\x )<0 and 
£(x*) = (3^/,(j;)/5/|^^^^. <0) (Figure 2D) (see Text S3. 3). Thus a 

monomorphic population should converge towards x and become 
stably maintained with no invasion of nearby mutants. In addition, 
similar to case (ii), a population that has evolved to x cannot be 
invaded by cheaters possessing the nuU mutation (because 
m(x ,0)<0; see Text S3.4). This scenario is called the "Intermediate 
evolution" case. 

The theoretical analysis described above is consistent with 
numerical simulations, in which an initial population of non-fixing 
rhizobia (x = 0) Mis into one of three evolutionary outcomes: case 
(i) "No evolution" when the cost (c) is relatively stronger than the 
benefit (b) (Figure 2A, crosses and 2E), case (ii) "Maximum 
evolution" when the benefit (A) is relatively stronger than the cost 
(c) (Figure 2A, squares and 2F), and case (iii) "Intermediate 
evolution" when the benefit and cost are balanced (Figure 2A, 
circles and 2G). 

3.1.4 Speed of the symbiosis evolution. In case (ii), an 
initial population of non-fixing rhizobia (x = 0) evolves to the 
maximum nitrogen fixation {x= 1) (Figure 2F). The evolutionary 
speed at which the symbiosis establishes is also affected by the 
cost-benefit balance. The lower the cost (or greater the benefit), 
the faster the evolution towards a cooperative behavior (Figure 3). 

3.1.5 Effect of the nonlinearity of the benefit 
function. Nonlinearity of the benefit function (ijv) does not 
qualitatively affect evolutionary behaviors, but influences the 
parameter conditions such that, as bj^ increases, the parameter 
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Figure 2. Effect of benefit and cost, assuming linear cost function. (A) Theoretically, the linear cost function (cn = 0) yields three evolutionary 
outcomes: (i) "No evolution" (gray), (ii) "Maximum evolution" (magenta), and (iii) "Intermediate evolution" (blue). This prediction is consistent with 
numerical simulations: cases (i), (ii) and (iii) correspond to crosses, squares and circles, respectively. (B-D) Pairvj/ise invasibility plots in (B) case (i), (C) 
case (ii), and (D) case (iii). Rare mutants with strategy y can invade the resident population with x in the gray region (i.e. w(x,y)>0), but cannot in the 
white region (i.e. w(x,y)<0). (E-G) Evolutionary dynamics of strategy distribution in (E) case (i), (F) case (ii), and (G) case (iii); darker shades indicate 
higher frequencies of a strategy. Parameters are: fa = 3.0 (B-G), fa/v = 0.0, c = 0.7 (B and E), 0.1 (C and F), and 0.22 (D and G), Cw = 0.0, n = 5. 
doi:1 0.1 371 /journal.pone.0093670.g002 



region of case (ii) reduces (Figure S 1 , magenta) while that of case 
(iii) expands (blue). This result signifies an inhibitory effect of bf^. 

3.2 Linear Benefit Function 

Next, we impose a linear benefit function, namely, hjf = 0. 
Similar to section 3.1, three cases of cooperative behavior emerge; 
(i) "No evolution" (Figures 4A and S2, gray), (ii) "Maximum 
evolution" (magenta), and (iii) "Intermediate evolution" (blue) (for 
details see Text S4). This model admits an intriguing outcome 
whereby nitrogen-fixing and cheating rhizobia coexist under the 
intermediate condition between cases (ii) and (iii) (orange, purple, 
and green). The remainder of this section is dedicated to this 
phenomenon. Coexistence emerges via two pathways: "evolution- 
ary branching" (orange and purple) and null mutation (green). 
Furthermore, the emergent cheating rhizobia are divided into two 
types depending on their relationship to the symbiosis system: "co- 
dependent" (orange) and "parasitic" (purple and green), in which 
nitrogen-fixing rhizobia (or the symbiotic system) are stabilized by 
"co-dependent" cheaters, but persist independently of "parasitic" 
cheaters. Thus, coexistence can be classified into the following 
three cases (for details see Text S4). 

3.2.1 Case (iv) "Co-dependent coexistence". In the 
orange parameter area in Figures 4A and S2, rhizobia adopt a 
singular strategy 0<x <1, which is CS but not ESS-stable (i.e. 
D'{x )<0 and E(x )>0) (Figure 4B). This situation is known to 



induce "evolutionary branching" or sympatric speciation, in 
which the monomorphic x population is invaded by nearby 
mutants and subsequently splits into two subpopulations with 
higher and lower activities of nitrogen fixation than the original x . 

This theoretical prediction is confirmed by numerical simula- 
tions, in which an initial population with x = 0 evolves to the 
singular strategy (x ) and then splits into two subpopulations 
(Figure 4A, open circles and 4E). These subpopulations gradually 
diverge until one attains the maximum (x= 1; fuU cooperator) 
while the other reduces to the minimum (x = 0; full cheater) under 
most of the investigated parameter conditions. Consequently, 
nitrogen-frxing and cheating bacteria evolve to stably coexist. 

This coexistence is predicted to be a co-dependent relationship, 
in which cooperators and cheaters require each other for their own 
survival because each monomorphic population is invaded by 
nearby mutants (i.e. Z)(0)>0 and D(1)<0; see Text S4). This 
prediction is also realized in numerical simulations. When 
cheating rhizobia are completely removed from the stable co- 
existence (Figure 4E, arrowhead), the remaining nitrogen-fixing 
bacteria lose their activities and their population destabihzes. This 
result suggests that cheating rhizobia stabilize the symbiotic 
relationship. Accordingly we name this type of coexistence "Co- 
dependent coexistence". 

3.2.2 Case (v) "Parasitic coexistence caused by 
evolutionary branching". In the purple parameter area in 
Figures 4A and S2, a singular strategy 0<x <1 exists that is CS 
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Cost (c) 

Figure 3. Evolution speed of symbiosis establisiiment. In case (ii) 
"Maximum evolution", an initial population of non-fixing rhizobia (x= 0) 
evolves to one containing full cooperators (x=1). The speed of this 
evolution increases (i.e. cycles required for "Maximum evolution" 
decrease), as the benefit (6) increases or as the cost (c) reduces. Each 
data point is the average of five calculations. Parameters are: 6a/ = 0.0, 
Ca, = 0.0, n = 5. 

doi:1 0.1 371/journal.pone.0093670.g003 

but not ESS-stable (Figure 4C), similar to ( asc Accordingly, 
cooperators and cheaters can co-evolve througli evolutionary 
branching. However, unlike case (iv), cooperators should persist 
regardless of the presence or absence of cheaters (because D{1)>0). 
This prediction is also confirmed by numerical simulations, in 
which cooperator bacteria can stably exist even after removal of 
the cheaters (Figure 4A, diamonds and 4F). In this case, the 
symbiotic relationship can persist independendy of the existence of 
cheating rhizobia. We refer to this type of coexistence as "Parasitic 
coexistence". 

3.2.3 Case (vi) "Parasitic coexistence caused by null 
mutation". In the green parameter region in Figures 4A and 
S2, the maximum strategy will evolve similarly to case (ii) (because 

D{x)>0 for 0<^<1). However, unlike case (ii), the evolved 
population of nitrogen-fixing symbionts {x—l) will likely be 
invaded by cheaters (j> = 0) possessing the nuU mutation (because 
Z£)(1,0)>0) (Figure 4D). 

This type of coexistence is also confirmed by numerical 
simulations. Nitrogen-fixing rhizobia first evolve as in case (ii), 
followed by transient introduction of a null mutation, which 
produces cheating rhizobia (Figure 4A, open squares and 4G). To 
investigate the fate of the cheaters existing in the population, we 
then remove this mutation. As theoretically predicted, cheating 
rhizobia can stably persist after ceasing the null mutation under 
the parameter conditions of case (vi). By contrast, under the 
conditions of case (ii), cheaters disappear immediately after 
mutation removal (Figure 4A, filled squares and 4H). 

3.2.4 Effect of the nonlinearity of the cost 
function. Nonlinearity in the cost function [cj^ promotes 
symbiosis evolution, because the parameter regions of case (ii) 
and (iii) increase and decrease, respectively, as c,v increases (see 
Figure S2, blue and magenta, respectively). Furthermore, Cjv 
induces the emergence of cheaters, because the parameter area in 
which cheaters can stably exist (cases (iv)-(vi)) is absent under weak 



nonlinearity of Cff<c^ = b/(n + nb + 2b), but extends with 
increasing Cj^ (Figure S2, orange, purple, and green). 

3.3 Nonlinear Benefit and Cost Functions 

Finally, we examine the effects of nonlinear benefit and cost 
functions. Similar to section 3.2, numerical calculation yields six 
evolutionary behaviors depending on the properties of the singular 
strategy. Under this generalized condition, the results are similar 
to those obtained in sections 3.1 and 3.2 (Cjv=0 and bjy=0, 
respectively). First, symbiotic evolution depends on the cost-to- 
benefit ratio in the following order: (i) "No evolution" (Figure 5, 
gray), (iii) "Intermediate evolution" (blue), (iv)-(vi) "Coexistence of 
cooperators and cheaters" (orange, purple, and green), and (ii) 
"Maximum evolution" (magenta). That is, the stronger the benefit 
(b) relative to the cost (c), the stronger the symbiotic relationship. 
Furthermore, Ajv inhibits the evolution of nitrogen-fixing bacteria; 
for example, the parameter region of case (ii) tends to decrease as 
ijv increases. In contrast, Cjv promotes the emergence of cheating 
rhizobia, as shown by the expansion of cases (iv)-(vi) as becomes 
stronger. 

3.4 Cost-Benefit Balance 

As discussed above, the cost-benefit balance determines the 
evolutionary' behavior of the symbiotic system. Here we summa- 
rize the evolutionary outcomes obtained in our model (Figure 6). 
According to the selection gradient D{x) (which defines the 
direction of evolution) (Figure 6, black arrows), an initially 
monomorphic population wiU evolve into three distinct nitrogen 
fixation activities: (A) full cheater with x = 0 if the cost is sufTiciently 
strong (case (i)), (B) full mutualism with x = 1 if the cost is weakened 
(cases (ii)/(vi)), and (C) partial mutualism with x — x at interme- 
diate cost (cases (iii)-(v)). The evolution of these monomorphic 
populations is affected and classified by their various evolutionary 
features. The emergence of a fully mutualistic population depends 
on whether cheaters evolving by null mutation can invade the 
population (^(1,0): invasibUity of cheaters). That is, a population of 
mutualists is stably maintained without invasion under weak cost 
conditions (!»(1,0)<0; case (ii)). As the cost is strengthened, 
mutualists and cheaters can coexist because cheaters can 
successfully invade under this condition (w(l,0)>0; case (vi)). 

However, the emergence of a partial mutualistic population 
depends on the occurrence of "evolutionary branching" (or 
sympatric speciation) (i.e. Efy*): ESS-stability of partial mutualists; 
see Figure 6, gray arrows). A partially mutualistic population 
persists without invasion by mutants under relatively strong cost 
conditions {E{x )<0; case (iii)). As the cost is weakened, evolution- 
ary branching creates a population of co-existent mutualists and 
cheaters (£(*:*)>0; cases (iv)/(v)). The cheaters in this type of 
coexistence can estabfish either a co-dependent or parasitic 
relationship with their mutualistic counterparts (i.e. i)(l): ESS- 
stability of mutualists). In the absence of cheaters, a population of 
mutualists can stably exist under weaker cost {D(\)>0; case (v)) but 
cannot maintain nitrogen fixation activity under stronger cost 
(Z)(1)<0; case (iv)). 

3.5 Efficiency of Nitrogen Fixation 

3.5.1 Assuming a linear cost function. The efficiency of 
nitrogen fixation (x^ relates to the average strategy adopted by the 
rhizobial population. Because a monomorphic population evolves 
in cases (i)-(iii), the efficiency in the linear cost function is 
described by 
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Figure 4. Effect of the benefit and cost, assuming a linear benefit function. (A) Theoretically, the benefit function (b/u = 0) yields six 
evolutionary outcomes: (i) "No evolution" (gray), (ii) "Maximum evolution" (magenta), (iii) "Intermediate evolution" (blue), (iv) "Co-dependent 
coexistence" (orange), (v) "Parasitic coexistence by evolutionary branching" (purple), and (vi) "Parasitic coexistence by null mutation" (green) (for 
derails see Text S4). This prediction is consistent with numerical simulations; crosses, squares, closed circles, open circles, diamonds, and open squares 
correspond to cases (i), (ii), (iii), (iv), (v) and (vi), respectively. (B-D) Pairwise invasibility plots of (B) case (iv), (C) case (v), and (D) case (vi). Rare mutants 
with strategy y can invade the resident population with x in the gray region (i.e. w(x,y)>0), but cannot in the white region (i.e. iv(x,y)<0). (E-H) 
Evolutionary dynamics of strategy distribution in (E) case (iv), (F) case (v), (G) case (vi), and (H) case (ii); darker shades indicate higher frequencies of a 
strategy. Once cheating bacteria with x<0.2 are removed from the coexistence situation (arrowheads), the remaining nitrogen-fixing bacteria can 
persist stably in case (v), but lose their activities in case (iv). A population of nitrogen-fixing rhizobia can be invaded by cheaters carrying the null 
mutation (brackets) in case (vi), but not in case (ii). Parameters: fa = 5.0 (B-H), c = 0.43 (B), 0.37 (C), 0.3 (D), 0.42 (E), 0.34 (F), 0.28 (G), and 0.26 (H). In all 
cases, 6^ = 0.0, C/v = 0.35, and n = 5. 
doi:1 0.1 371 /journal.pone.0093670.g004 



1 case(ii) (c<ci) 
Xeff=<. X* case(iii) (ci <c<co) 
0 case(i) (cq < c) 



(8) 



where x is described in Text S3. 3. The efficiency continuously 
decreases as the cost c increases (Figure 7 A). 

3.5.2 Assuming a linear benefit function. When full 
cooperators {x=l) and full cheaters (* = 0) coexist in cases (iv)- 
(vi), the efficiency of nitrogen fixation in the linear benefit function 
depends on the proportion of cooperators, i.e. 

p* = {h-(n + h){\-CM)c}/in-m- CM)bc (9) 
(for details see Text S4.8). The efficiency then becomes 



1 case(ii) 

p* cases(iv) (vi) 

X* case(iii) 

0 case(i) 



(10) 



where x is described in Text S4.4. As observed in the linear cost 



function, reducing the cost c increases the efficiency, but an 
interesting discontinuity occurs at the transition between case (i)/ 
(iii) and case (ii)/(iv)/(v)/(vi) (Figure 7B-D). 

3.6 Effect of Root Nodule Number 

Overcrowding of nodules on a host is expected to impede the 
mutualistic relationship, because the effect of an individual 
rhizobium on the benefit function B(x) weakens as the nodule 
number (n) increases (Eq. (1)). Supporting this prediction, the 
parameter conditions of case (i) expand with increasing n (Figure 8, 
gray) while those of cases (ii) and (iii) reduce (magenta and blue). 
Furthermore, higher nodule numbers promote the emergence of 
cheaters; the parameter area occupied by cases (iv)-(vi) is absent at 
n = 1 but expands as n increases (orange, purple, and green). These 
results suggest that, as nodule number increases, the symbiotic 
relationship destabilizes and is less easily established, while 
cheating rhizobia readily emerge. 

3.7 Effect of Partner Choice 

Partner choice is often described as the preferential selection of 
colonizers by the host plant [9-11]. To examine the effect of 
partner choice, we introduce a parameter oc(x): the probability of 
accepting rhizobial colonization (where OSo:(x)<l and oc'[x)>0). 
We then consider the invasibility of rare mutants with nitrogen 
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0.0 ^ 5.0 0.0 ^ 5.0 0.0 ^ 5.0 



0.0 0.1 0.2 0.3 0.4 
Nonlinearity of benefit 



0.5 



Figure 5. Effects of nonlinear benefit and cost functions. Nonlinear benefit and cost functions yield the evolutionary outcomes (i) "No 
evolution" (gray), (ii) "iVIaximum evolution" (magenta), (iii) "Intermediate evolution" (blue), (iv) "Co-dependent coexistence" (orange), (v) "Parasitic 
coexistence by evolutionary branching" (purple), and (vi) "Parasitic coexistence by null mutation" (green). Parameter regions of these cases are 
delineated by the properties of their singular strategy (similar to Text S4); case (i) D(0)<0, case (ii) D(x)>0 for 0<x<1 and w(1,0)<0, case (vi) D(x)>0 
for 0<x<l and w(1,0)>0, case (iii) x* is CS and ESS-stable (i.e. D'(x*)<0 and f(x*)<0), case (iv) x* is CS but not ESS-stable (i.e. D'(x*)<0 and f(x")>0) and 
an unstable monomorphic population of cooperators exists (i.e. D(1)<0), and case (v) if x is CS but not ESS-stable (i.e. D'(x )<0 and f(x )>0) and a 
stable monomorphic population of cooperators exists (i.e. D(1)>0), where x is the smallest singular strategy (D(x ) = 0 and 0<x <1). In all cases, n = 5. 
doi:1 0.1 371/journal.pone.0093670.g005 



fixing ability y into the resident population with x. The mutant 
colonizes the host at a a()j)/o£(.v-)-fold different rate than in the 
absence of partner choice (Figure IB). Equivalently, partner choice 
affects the fitness of the mutant (^'xCi')) by a factor of aiy)! a{x), so 
that 



Fx(y)=fxiy)a.(y)/a(x). 

Accordingly, we obtain: 

D{x) = dFx(y)/dy\^^^ = D,(x)/a{x), 

D'{x*) = D'^{x*)/a{x'), 



E{x) = 8^F,(y)/dy^ 



--E^{x)/ix(x), 



(11) 



(12a) 



(12b) 



(12c) 



where D,(x) = dfjdyl^^, E,(x) = d^f,/dy^\^,^^, 

fx=fx(y)cx.(y) = B{x)C(y)ix(y), and x* is an evolutionarily singular 
strategy (i.e. D{x ) = 0). This result indicates that a singular strategy 



(satisfying D(x ) = 0) and its CS and ESS stabilities (which depend 
on the signs of Z)'(x) and E{x), respectively) can result from a 
modified fitness function ^, = B(x)CJXi), where Ca(x) = C{x)a(x) 
is a modified cost function. Therefore, similarly to host sanction, 
partner choice promotes symbiotic evolution by inffuencing the 
cost function. 

3.8 Effect of Mixed Nodule Populations 

Mixed nodules colonized by multiple symbionts are common in 
the laboratory [40-42] and are also observed under field 
conditions [43]. Mixed nodule populations are thought to reduce 
the evolutionary effects of host sanction and thereby promote the 
emergence of ineffective rhrzobia [10,13,16,34]. To examine the 
effects of mixed nodule populations in our model, we introduce a 
parameter p (0£/?£l), specifying the frequency of mixed nodules 
(or doubly colonized nodules). Each nodule is colonized by two 
rhizobial strains with probability ;8; otherwise it is colonized by a 
single strain (with probability 1-/J). 

In a mixed nodule, two strains adopting strategies xi and X2 
would proliferate at rates depending on their relative fitness C{xi) 
{i= 1 or 2). Thus, we assume their proportion in the population as 
Pi = C(Xi)/(C(xi) + C(x2)). Accordingly, the nitrogen frxation 
activity of the mixed nodule, which contributes to the average 
nitrogen frxation (x) in the benefit function, is given by 
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(i) 
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(iii) 

Intermediate 
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0 A- 1 

Nitrogen fixation 

(X) 

(iv) 

Co-dependent 
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Evolutionary 
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D(.r) > 0 
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Parasitic 
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Figure 6. Cost-benefit balance in tlie symbiosis evolution. Cost-benefit balance determines the evolution of the symbiotic system by 
affecting various evolutionary features, such as selection gradient D(x), w(1,0): invasibility of cheaters in a population of mutualists, E(x'): ESS-stability 
of partial mutualists, and D(1): ESS-stability of mutualists. Thereby, the evolutionary outcomes obtained in our model can be classified into six cases 
(i)-(vi) according to the cost-benefit balance. For details see text. The selection gradient determines the direction of evolution, such that a 
monomorphic population evolves towards larger strategies if D(x)>0 but towards smaller strategies if D(x)<0 (black arrows). Circles indicate singular 
strategies that are CS (i.e. D(x)=0 and D'(x)<0). Filled squares correspond to cheaters (x = 0) or mutualists (x=1) that are locally ESS-stable (i.e. 
D(0)<0 or D(1)>0, respectively), and open squares correspond to those that their strategy is not ESS-stable but can coexist with the other strategy. 
doi:1 0.1 371 /journal.pone.0093670.g006 



piXi +P2X2. Furthermore, the fitness of each strain depends on its 
proportion (p-j in tlie mixed nodule and is simply defined as 
B{x)C(xi)pi, which is reduced to single colonization (i.e. 
B(x)C(Xi)) and 0 when p, = 1 and p, = 0, respectively. 

Now consider that mutants with nitrogen fixation ability ji 
invade a resident population with ability x. If the mutant is rare, 
the probability that the host plant is colonized by multiple mutants 
is very low and can be neglected. Thus, root nodules will harbor 



single y colonies and colonies containing x and y with relative 
frequencies of (1-/J) and 2/?, respectively (Figure IC). The fitness of 
the mutant in a single colonization is fi=B{x)C(j') where 
x = {{n—\)x-\-y)/n, and that in a double colonization is 
f2 = B(x)C(y)py, where x = {(n—l)x + z)/n, z=p^-x+pyy, 
p, = C(x)/{C{x) + C(y)), and Py = l -p, = C(y)/(C(x) + C(y')). 
Consequendy, the expected fitness of mutant j is 




Figure 7. Efficiency of nitrogen fixation. The efficiency of nitrogen fixation (x^^) decreases with increasing cost c. (A) This decrease is continuous 
for a linear cost function (Cfj = 0). (B-D) If the benefit function is also linear (b^, = 0), the decrease is discontinuous (arrowheads) at the transition 
between cases (i) and (ii) (B), cases (i) and (v) (C), and cases (iii) and (v) (D). The parameter regions of cases (i), (ii), (iii), (v) and (vi) are indicated in gray, 
magenta, blue, purple, and green, respectively. Parameters are: b = 3.5 (A), 1.0 (B), 2.0 (C), and 4.0 (D); b„ = 0.2 (A) and 0.0 (B-D); c„ = 0.0 (A) and 0.5 (B- 
D); n = 5. 

doi:1 0.1 371 /journal.pone.0093670.g007 
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Figure 8. Effect of root nodule number. As the nodule number (n) 
on a host root increases, the symbiotic relationship evolves less easily. 
The parameter regions of case (ii) (magenta) and case (ill) (blue) 
decrease while that of case (i) (gray) increases. However, cheating 
rhizobia emerge more easily, as shown by the expanding parameter 
region in which cooperators and cheaters coexist (cases (iv)-(vi); 
orange, purple, and green). Theoretical and numerical predictions are 
indicated respectively by gray area and crosses (case (i)), magenta area 
and closed squares, (case (ii)), blue area and closed circles (case (ill)), 
orange area and open circles (case (iv)), purple area and diamonds (case 
(v)), and green area and open squares (case (vi)). Parameters are: b = 5.0, 
b =0.0, Cjv = 0.4. 
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F,(y) = ((l-p)fi+2ftf2)/(i+f^)- 



Thus we obtain the selection gradient 



(13) 



--Do(x)(\- p/2)/{l+ p) + B{x)C'(x)p/(l + 1^1 



(14) 



where Di}(x) = dfi(y)/dy\^,^^ is the gradient in a pure population 
(i.e. P = 0). The second term on the right-hand side of Eq. (14) is 
always negative because C'(x)<0, implying that mixed nodule 
populations inhibit the evolution of nitrogen frKation. 

This prediction is supported by numerical simulations, in which 
the parameter region of "No evolution" expands as P increases, 
while that of "Maximum evolution" decreases (Figure 9, gray and 
magenta). In contrast, mixed nodule populations exert no positive 
effects on mutualistic-cheater co-existence, because the parameter 
range of cases (iv)-(vi) is not largely influenced by P (Figure 9, 
orange, purple, and green). However, as P increases, this 
parameter region moves toward lower cost (c). In this scenario, 
cheating rhizobia will either emerge or disappear in mixed nodule 



populations, depending on the parameter conditions. For example, 
setting c = 0.2 in Figure 9, cheaters cannot exist in mixed nodule 
populations that establish with low probability (/?</?] —0.2 17; case 
(ii), magenta), but can emerge when the probability of mixed 
populations is higher (/J>/?i; cases (iv)-(vi), green, purple, and 
orange). Conversely, when c = 0.4, cheaters can coexist with 
mutualists at lower P (</?2~0.212; cases (iv) and (v), orange and 
purple), but are excluded at higher P {>p2', cases (i) and (iii), blue 
and gray). 

Discussion 

In .symbiotic relationships, the participating organisms provide 
mutual benefits to each other. Such positive fitness feedback 
reinforces their mutualistic interaction. However, symbiotic 
systems encourage the emergence of selfish parasitic cheaters, 
whose performance undermines the system. Thus, symbiotic 
systems are simultaneously exposed to promoting and destabilizing 
forces, analogous to benefit and cost in game theory. 

One of the most famous symb loses occurs between legumes and 
rhizobia. Here, we intensively investigated how benefit and cost 
influence the evolution of this symbiosis, and the conditions 
required for establishing the symbiotic relationship. According to 
our model, stable mutualism depends on the cost-benefit balance 
(Figures 6 and 10). That is, a tight symbiotic relationship emerges 
when the beneficial effect is much stronger than the cost (case (ii)). 




0.0 ' 



D(0) = 0 
£(.v-) = 0 
D(1) = 0 
clet = 0 
if(l,0) = 0 



0.0 0.2 0.4 0.6 0.8 1.0 

Frequency of mixed nodule (J^) 

Figure 9. Effect of mixed nodule populations. Our model predicts 
that mixed nodule populations destabilize the symbiotic relationship 
(see section 3.8 for details). This prediction is supported by numerical 
simulations. Theoretical and numerical predictions are indicated 
respectively by gray area and crosses (case (i)), magenta area and 
closed squares, (case (ii)), blue area and closed circles (case (iii)), orange 
area and open circles (case (iv)), purple area and diamonds (case (v)), 
and green area and open squares (case (vi)). Parameters are: b = 4.0, 
faw = 0.0, Cn = 0.35, n = 5. 
doi:10.1371/journal.pone.0093670.g009 
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but is dissolved under the opposite condition of relatively strong 
cost (case (i)). In the intermediate condition, where benefit is 
approximately offset by cost, more complicated behaviors emerge 
such as imperfect symbiotic interactions (case (iii)) and the 
coexistence of cooperators and cheaters (cases (iv)-(vi)). As the 
benefit strengthens relative to the cost, the evolutionary outcome 
shifts in the order of "No evolution", "Intermediate evolution", 
"Coexistence of cooperators and cheaters", and "Maximum 
evolution" (Figure 10). 

Ineffective rhizobia with little or no nitrogen frxation activity are 
widely distributed in natural environments [5-8], and their 
emergence and persistence has aroused much discussion 
[4,10,44^6]. Our model showed that cheaters can emerge when 
the benefit and cost are balanced (Figure 10, cases (iv)-(vi)). Under 
such conditions, nitrogen-fixing cooperators can co-exist with 
cheaters and maintain stable symbiotic relationships. Therefore, 
the emergence of cheating bacteria does not necessarily disrupt a 
symbiotic relationship. 

According to our model, cheating rhizobia can be classified into 
"co-dependent" or "parasitic", which exert different effects on the 
system. The presence of "co-dependent" cheaters stabilizes a 
symbiotic system (Figure 4E), while "parasitic" cheaters produce 
no effect on the system (Figure 4F, G). However, in natural 
environments, these classes of cheaters may be difficult to 
distinguish. One reason is that cheaters are spontaneously 
generated by mutations; thus, the mutualistic effects of cheaters 
in the complete absence of the nuU mutation are not easily 
examinable. In addition, cheaters could easily transfer between 
"co-dependent" and "parasitic" under small changes of parameter 
conditions, because the classes are adjacent in the parameter space 
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Figure 10. Model for the evolution of the legume-rhizobia 
symbiosis. The evolution of the legume-rhizobia symbiosis depends 
on the cost-benefit balance. As the benefit strengthens relative to the 
cost, the evolutionary outcome shifts in the following order: (i) "No 
evolution", (iii) "Intermediate evolution", (iv)-(vi) "Coexistence of 
nitrogen-fixing and cheating rhizobia", and (ii) "IVlaximum evolution". 
The symbiotic relationship is reinforced by partner fidelity feedback, 
which strengths the benefit, and by host sanction and partner choice, 
which diminish the cost. In addition, as the number of nodules on a 
root increases, symbiotic rhizobia are displaced by selfish cheaters. 
doi:10.1371/journal.pone.0093670.g010 



(Figures 4A, 8, 9, and S2). Because parameter values can be 
perturbed by environmental fluctuations, cheating rhizobia may 
continuously alternate between "co-dependent" and "parasitic". 

As mentioned above, the persistence of the legume-rhizobia 
symbiosis in nature, despite the ubiquity of ineffective rhizobia, has 
aroused much interest. Various mechanisms that stabilize the 
symbiotic relationship have been proposed; in particular, partner 
fidelity feedback [1-4], partner choice [4,9-12], and host sanction 
[4,10,13-16]. In our model, partner fidelity feedback provides a 
benefit (or promoting force), while host sanction and partner 
choice reduce the cost function C(x) (or destabilizing force). 
Therefore, our model suggests that these factors strengthen the 
symbiotic interaction in opposite ways; partner fidelity feedback 
reinforces the benefit while host sanction and partner choice 
ameliorate the cost (Figure 10). Thus, although the legume- 
rhizobia symbiosis may be maintained by a single mechanism, it is 
rather more likely to be cooperatively reinforced by various 
stabilizing mechanisms. 

Mixed nodule populations are thought to reduce the evolution- 
ary effects of host sanction, and thereby encourage the persistence 
of cheating rhizobia [10,13,16,34]. However, our model predicts 
that mixed nodule populations exert destabilizing effects on the 
symbiotic relationship, and can cause cheaters to either emerge or 
vanish depending on the parameter conditions. This evidence 
suggests that mixed nodule populations produce both positive and 
negative effects on the emergence of cheating rhizobia in different 
situations. 

In this paper, we assumed a decreasing cost function (i.e. 
C(x)<0), implying that more cheating rhizobia proliferate more 
efificiendy if they infect the same host plant, because the fitness is 
given by =5(x)C(x,) (Eq. (1)). Interestingly, under these 
conditions, ineffective rhizobia do not necessarily evolve to 
persistence, and tight mutualistic interaction is not prohibited 
(for example case (ii)). Thus, while cheaters are benefitted by a 
growth advantage, they do not necessarily co-evolve with 
cooperators or cause collapse of the symbiotic system. This result 
reconfirms that the cost-benefit balance, rather than the cost 
alone, is important for symbiotic evolution (Figure 10). 

However, double inoculation experiments have shown that 
nitrogen-fixing strains proliferate with equal or superior efficacy to 
ineffective strains colonizing the same plant seedling 
[4,11,15,33,47-50]. This experimental observation may be 
attributed to various synergistic effects, but ultimately suggests 
that the cost function is not a simple decreasing function of 
nitrogen fixation. Under such conditions (i.e. C'(x)SO), our theory 
predicts that rhizobial populations always evolve to the maximum 
activity of nitrogen fixation (see Text S2). However, this prediction 
contradicts the ecological fact that ineffective rhizobia are 
widespread in natural environments [5-8]. Possibly, the benefit 
and cost functions are not simple increasing or decreasing function 
of X as we have assumed, but respond with more complexity to 
nitrogen resources. Furthermore, our highly-simplified model 
ignores a number of effects such as plant strategy dynamics, 
symbiont growth in the soil, spatial heterogeneity, and environ- 
mental fluctuations. Incorporation of these effects into our model 
might introduce greater complexity to the evolution of legume- 
rhizobia symbiosis. 

In our model, rhizobia fitness is influenced by both the average 
nitrogen fixation of the rhizobia infecting a host plant {x) and 
individual rhizobium activity (x,) (Eq. (1)), here denoted by 
"systemic" and "local" effects, respectively. Both effects are 
essential for coexistence of nitrogen-fixing and cheating rhizobia, 
because coexistence cannot emerge without their mutual cooper- 
ation (for details see Text S5). Although the benefit and cost 



PLOS ONE I www.plosone.org 



11 



April 2014 I Volume 9 | Issue 4 | e93670 



Evolutionary Dynamics in Legume-Rhizobia Symbiosis 



functions in this paper are simply given by the systemic and local 
eflects, respectively, they could take more complex general forms 
incorporating both effects; for instance, B(x,Xi) and C(x,Xi). If the 
systemic and local effects in these functions are approximately 
separable from each other (i.e. B(x,Xi)^ Bs(x)Bi^(Xi) and 
C{x,Xi) ~ Cs(x)Cl(x,)), the fitness can be written as 

(Pi = B(x,xi)C(x,Xi) ~ S{x)L(x,), (15) 

where S{x) = Bs(x)Cs(x) and L(x) = Bl(x)Cl(x) are functions of 
the systemic and local effects, respectively. Replacing S(x) and L{x) 
respectively with B(x) and C{x), Eq. (15) takes the fitness form of Eq. 
(1); hence, the evolutionary dynamics described by Eq. (15) can be 
understood within the framework described in this paper. For 
example, the benefit (or partner fidelity feedback) might include 
both local and systemic effects. The local beneficial effect (i.e. 
-8z,(^t)) would eventually reinforce the symbiotic relationship by 
affecting the local function L{x} of Eq. (15), similarly to the 
weakening of the cost function Qx] by host sanction and partner 
choice in Eq. (1). 

In this paper, we constructed and investigated the evolutionary 
dynamics of the legume-rhizobia symbiotic system, driven by the 
opposite effects of promoting force (or benefit) and destabilizing 
force (or cost). We then determined a comprehensive set of 
conditions under which the symbiotic relationship is established 
and cheating bacteria emerge. Our findings provide a theoretical 
basis for understanding how the legume /rhizobium symbiosis 
evolves. However, as mentioned above, our model is extremely 
simpKfied, and evolutionary behaviors in the natural environment 
are expected to be much more complicated than those seen in our 
model. By incorporating such effects into our model, we could 
better understand the evolution of symbiotic systems. 

Supporting Information 

Figure SI Effect of b^, assuming a linear cost function. 

If a linear cost function (cjtf= 0) is assumed, nonlinearity in the 
benefit function (bj^ hinders symbiotic evolution. As increases, 
the parameter region of case (ii) "Maximum evolution" (magenta) 



decreases while that of case (iii) "Intermediate evolution" (blue) 
increases. This prediction is consistent with numerical simulations; 
cases (i), (ii) and (iii) are indicated by crosses, squares, and circles, 
respectively. Parameters are: b = 5.0, cj^= 0.0, n = 5. 
(TIF) 

Figure S2 Efiect of Cj^, assuming a linear benefit 
function. If a linear benefit function (A v = 0) is assumed, 
nonlinearity in the cost function (cjy) promotes symbiosis evolution. 
As Cjv increases, the parameter region of case (ii) "Maximum 
evolution" (magenta) increases, while that of case (iii) "Interme- 
diate evolution" (blue) decreases. In addition, Cjif promotes the 
emergence of cheaters. The parameter region in which nitrogcn- 
frxing and cheating rhizobia coexist (cases (iv)-(vi); orange, purple, 
and green areas, respectively) is absent when «jv<Ov - b/(ri+(n+2) 
b), but otherwise increases as cjv increases. These predictions are 
consistent with numerical simulations; cases (i), (ii), (iii), (iv), (v) and 
(vi) are indicated by crosses, closed squares, closed circles, open 
circles, diamonds, and open squares, respectively. Parameters are: 
A = 5.0, Ajv=0.0, re = 5. 
(TIF) 

Text SI Theoretical framework of adaptive dynamics. 

(PDF) 

Text S2 Increasing cost function (C(«)>0). 

(PDF) 

Text S3 Linear cost function (cjv= 0). 

(PDF) 

Text S4 Linear benefit function (&iv= 0). 

(PDF) 

Text S5 Systemic and local effects of nitrogen fixation. 

(PDF) 
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